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ABSTRACT 


An  analytical/numerical  procedure  has  been 
developed  which  yields  the  optimum  geometry  for  a 
constrained  hole  shape  in  a  large  plate  under 
prescribed  boundary  stresses  at  infinity.  The 
optimality  criterion  is  based  on  the  minimization 
of  a  certain  stress  integral  taken  around  the  hole 
boundary.  Muskhelishvili' s  method  is  used  to 
first  obtain  the  symbolic  stress  expressions  for  a 
given  mapping  function  with  unknown  coefficients. 
These  stress  expressions  are  then  squared  and 
integrated  around  the  hole  boundary  to  obtain  the 
stress  integral.  The  squaring  prevents  any  cancel¬ 
lation  of  compressive  and  tensile  stresses  during 
the  integration  process  and  thus  produces  a  more 
meaningful  integral. 

A  sample  problem  is  worked  out  in  detail  to 
demonstrate  this  procedure.  The  optimum  shape  of  a 
square-like  (double  barrel  shape)  hole  with  rounded 
comers  was  determined  by  restricting  the  mapping 
function  to  only  one  specific  unknown  coefficient. 
Modification  of  this  shape  by  introduction  of  an 
additional  term  is  also  discussed.  Numerical 
stress  concentration  values  for  this  case  are 
compared  with  those  from  other  sources. 

Most  of  the  algebraic  and  numerical  calcula¬ 
tions  presented  in  this  paper  were  performed  using 
MACSYMA,  a  symbolic  manipulation  package  developed 
at  MIT  and  in  regular  use  at  DTNSRDC. 


ADMINISTRATIVE  INFORMATION 

This  work  was  sponsored  by  the  Naval  Sea  Systems  Command  and  funded 
during  Fiscal  Year  1979  under  the  Mathematical  Sciences  Program,  Program 
Element  61153N,  Task  Area  SR0140301,  Work  Unit  1-1808-010. 

INTRODUCTION 

One  of  the  most  important  problems  in  the  design  of  plate  structures 
is  to  determine  and  minimize  the  stress  concentration  due  to  the 
inevitable  presence  of  holes  and  other  discontinuities.  The  literature 
contains  a  great  deal  of  information  on  this  age-old  problem;  however, 
most  of  this  information  is  concerned  with  direct  determination  of 
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stresses  for  given  geometry  and  load  conditions.  Some  recent  exceptions 

1* 

are  the  experimental  results  reported  by  Durelli  and  Rajaiah,  the  harmonic 
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hole  concept  discussed  by  Bjorkman  and  Richards,  ’  and  the  optimization 

4 

approach  developed  by  Schnack.  The  experimental  results  of  Durelli  and 
Rajaiah  were  obtained  by  step-by-step  machining  of  hole  boundaries  in  a 
photoelastic  model  until  the  stress  (related  directly  to  the  order  of  iso- 
chromatic  fringe  for  load  free  boundaries)  was  observed  to  be  approximately 
constant  in  both  the  tensile  and  compressive  boundary  regions.  This  tech¬ 
nique  is  effective,  but  it  does  call  for  a  series  of  experiments  for  each 
unique  situation.  The  harmonic  hole  concept  of  Bjorkman  and  Richard  con¬ 
serves  the  first  stress  invariant;  however,  it  appears  to  have  limited 
application  as  an  optimization  procedure  because  it  results  in  somewhat 
unpractical  hole  shapes.  Schnack 's  procedure,  an  interesting  finite 
element  approach  based  on  the  fade-away  law  developed  by  Neuber,  is  in  a 
way  the  numerical  equivalent  of  Durelli  and  Rajaiah 's  experimental  method. 

The  present  paper  is  a  step  toward  practical  design  optimization,  in 
that  the  technique  permits  the  designer  to  have  some  control  over  the  gen¬ 
eral  shape  of  the  hole.  The  problem  can  thus  be  formulated:  Determine  the 
geometry  of  a  given  type  of  hole  (e.g.,  square-like)  which,  when  placed  in 
a  uniformly  uni-  or  biaxially  loaded  large  plate,  produces  boundary  stresses 
which  minimize  a  certain  stress  integral.  Such  an  optimum  hole  appears  to 
lead,  in  most  practical  cases,  to  a  very  desirable  overall  stress  field. 

The  approach  involves  first  determining  the  symbolic  boundary  stress 
expression  for  a  desired  mapping  function  in  a  fairly  general  form 
following  the  well  established  Muskhelishvili  method.  This  step  may 
require  considerable  algebra  depending  on  the  mapping  function.  In  the 
present  case  most  of  these  algebraic  operations  and  the  following  steps 
were  accomplished  by  MACSYMA,  a  symbolic  manipulation  package  developed 
at  MIT.  The  next  step  is  to  square  this  expression  and  integrate  it 
around  the  hole  so  that  the  tensile  and  compressive  stresses  do  not 
annihilate  each  other  during  integration.  The  last  step  is  to  minimize 
this  expression  with  respect  to  the  mapping  function  coefficients.  These 


*  A  complete  listi  'g  of  references  is  given  on  page  13. 
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extreme  values  of  coefficients  contain  the  necessary  information  to 
define  the  geometry  of  the  desired  opening. 


STRESS  DETERMINATION 


Let  the  mapping  function  be  represented  by 


m  m,  m_  m 

z(0  =  s+^  +  ^  +  -|  +  -4 

c  c3  c5  c7 


a  function  which  can  map  a  unit  circle  in  ?-plane  (?  =  e  p)  into  a  hole 

with  at  least  two  areas  of  symmetry  in  the  z-plane  (z  =  x  +  iy).  The 

coefficients  m  (n  =  l,  3,  5,  7)  are  assumed  to  be  real  and  have  the  usual 
n 

limitations  on  them  which  assure  that  z'(?)  #  0  for  ?  *  1  (prime 
indicates  differentiation  with  respect  to  the  argument).  The  stress 
field  can  be  represented  by 


where 


o  -o  +  2ix  . 
8  a  a8 


P2*'(C) 


[zU)  *’(c)  +  ♦’(01 


<K?)  =  5  o  >  o„  and  x  „  are  the  usual  plane  stresses; 

z  (?)  a  8  aB 


p  is  the  absolute  value  of  ?;  and  bar  (-)  over  a  function  represents  its 
complex  conjugate.  The  functions  <)>(£)  and  i|»(?)  are  known  to  have  the 


♦  (c)  =  s  ?  +  i 


♦(C)  =  D  ?  +  E 


where  a  and  b  are  constants  to  be  determined  and 
n  n 


3 


s 

D 


(6) 


where  p  and  q  represent  the  uniform  stresses  at  infinity  at  an  angle  0 

to  the  x-  and  y-axes,  respectively.  The  a  were  determined  using 

n 

Muskhelishvili' s  procedure  for  a  load-free  hole  boundary  and  0=0,  and  are 
given  in  the  appendix  for  completeness  sake. 


OPTIMIZATION  PROCEDURE 


Since  the  hole  boundary  (p=l)  is  assumed  to  be  load-free 

(a^  =  t  =  0),  the  stresses  around  its  boundary  can  be  described  by 

equation  (2)  which  requires  the  determination  of  <Hc),  i.e.,  a  only. 

n 

Now  the  mapping  function  z(£),  given  by  Equation  (1),  can  be  redefined  in 

terms  of  Q  ,  such  that 
n 


m 


1 


\ 


m3  3  S 

i=l  j=i+l 


r  qiQJ 


c  2  2  E  Q.Q.Q. 

J  i=l  j=i+l  k=j+l  1  J  K 


m7  =  - 


Q2  Q3  Q4 


(7) 


These  substitutions,  which  can  be  reduced  to  a  cubic  equation  for 
determining  Q  ,  permit  z'(?)  to  be  given  by 

4  Qn 

z’(0  =  n  (l  -  )  (8) 

n=l  52 

Clearly  | j  <  1 ,  and  the  stresses  then  can  be  represented  at  the 
hole  boundary  by 


4 


(9) 


4  (2n-l)a 


4  Re 


S  -  I 
n=l 


2n-l 


2n 


n  (i  -  -r  ) 


n=l 


which  can  conveniently  be  simplified  to 


og  =  4  Re 


P+ 

L  n=1  (i--^)J 


(10) 


where  use  has  been  made  of  the  last  of  Equations  (7)  as  well  as  that  of 


the  value  of  a^  from  the  appendix.  The  P^  are  defined  as  follows: 


Pl  =  - 


Q1S  ~ Qlal  "  3Qla3  "  5Qia3  ~ 7a7 


q1(q2-q1)(q3-q1)(q4-q1) 


P2’  P3’  and  P4  Can  be  obtained  fay  cyclically  rotating  in  the  above 
expression.  Equation  (10)  reduces  to 


=  4[-S+  Z  (P1Q""1+P2Q2  1+p3Q3"1+p4QJ"1)cos  2(n-l)3]  (11) 

n=l 

Although  the  constants  Pn  and  can  be  individually  complex,  the 
expression  to  be  summed  in  Equation  (11)  will  always  be  real  for 
symmetrically  located  holes  and  0=0.  Parseval's  theorem  can  now  be  used 
to  show  that 

.  2ir  4  2  4  2  4  4  P  P 

1  "  dB  =  2(s  "  s  pn)  -V  Pn)  +  1  1  T~n~n  (12) 

1611  0  6  1  n  1  n  11  1-QmQn 

where  I  is  some  form  of  a  measure  of  the  strain  energy  around  the  hole. 

It  should  be  noted  that  the  integration  in  Equation  (12)  is  performed 
around  3  and  not  around  the  polar  angle  0 .  Equation  (12)  can  be 
differentiated  with  respect  to  and  equated  to  zero  to  minimize  1.  By 
introducing,  for  example,  some  form  of  Lagrange's  multipliers,  it  is 
possible  to  introduce  additional  constraints. 
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Whenever  more  than  one  constant,  m  or  Q  ,  is  involved  in  the 

n  n 

mapping  function,  the  expression  for  I,  Equation  (12),  explicit  in 
only,  becomes  so  complex  that  its  derivatives  with  respect  to  Qn  produce 
unwieldy  expressions.  For  this  reason,  it  may  be  desirable  in  such 
cases  to  simply  compute  the  numerical  values  of  I  over  a  range  of 
parameters  and  thus  determine  the  minimum  point . 

Equations  (7)  thru  (12)  have  been  developed  for  the  mapping  function 
given  by  Equation  (1).  Therefore,  these  expressions  should  be  directly 
reducible  to  accommodate  other  mapping  functions  with  any  of  the 
equal  to  zero.  While  this  is  generally  true,  some  expressions,  such  as 
the  cyclic  definitions  of  P  ,  appear  to  have  poles  in  some  situations. 
However,  these  are  not  real  poles  and  are  easily  removable  as  the  reader 
can  verify.  Special  attention  should  be  paid  when  programming  these 
general  expressions  for  a  computer. 


EXAMPLE  I 

Let  the  mapping  function  be  given  by 

m 

z  =  i,  +  —  ,  |m3|  <  -  ;  (13) 

which  for  non-zero  values  of  m^  can  map  a  unit  circle  into  a  square  with 
rounded  corners.  The  actual  geometry  of  the  opening  (e.g.,  rounded 
corners)  depends  on  the  value  of  m^.  Equations  (7)  and  (10)  can  be  used 
to  obtain 


Ql  =  -Q2  =  /3^ 


Q3  "  Q4  *  ° 


Px  =  S  + 


P2  =  S  + 


p  ■  p 
3  4 


2(3-Qj)Q1 


2(3-Q2)Q2 


The  expressions  for  Q  and  P  given  in  Equations  (14)  can  be  substituted 

n  n 

in  Equation  (11)  and  reduced  to 


°6  =  4 


(l-m^)  (l-9m3)  S  +  (l-3m^)D  cos  23 
(l-m^)  (l+9m2  -  6m3  cos  46) 


These  values  of  and  P  ,  when  substituted  in  Equation  (12),  yield  the 
integral  I  as 


I  =  -2S2  +  + 


l-9m2  (l-m3)2(l-9m2) 


The  first  derivative  of  this  integral  with  respect  to  yields  the 
following  condition  for  extremum  values  cf  I: 

36  S2  m^  -  108  S2m2  +  18(D2+6S2)m2  -  9(D2+4S2)m3~D2  =  0  (17) 

The  symbolic  solution  of  Equation  (17)  was  obtained  as  a  function  1 
S  and  D,  but  it  is  not  included  here  because  of  space  limitations.  Four 
specific  cases  were  considered: 

(i)  isotropic  loading,  p  *  q  *  1; 

(it)  uniaxial  loading,  p  =  1,  q  =  0; 

(iii)  biaxial  loading,  p  =  2,  q  =  1;  and 

(iv)  shear  loading,  p  =  1,  q  =  -1. 

The  values  of  m  which  were  found,  using  the  solution  of  Equation  (17),  to 
be  within  the  bound  of  Equation  (13)  (i.e.,  | j  <  —  for  the  four 
cases)  are,  respectively: 

(i)  m3  =  0; 

(ii)  m.  =  -.049297  (^  -0.05); 


(ii)  m3  =  -.049297  (^  -0.05); 

(iii)  m3  =  -0.0107713  (%  -0.01);  and 

(iv)  m3  =  -0.0936  (^  -0.09)  . 

Case  (i) 

The  first  case  is  obvious,  since  ra3  =  0  represents  a  circular 
opening  which  is  optimum  in  an  isotropic  stress  field  characterized  by 
p  «  q  =  1.  The  value  of  I  is  0.5. 


The  second  case  (m3  =  -0.05,  p  =  1,  q  =  0)  is  worthy  of  comment. 
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The  value  of  m^  used  by  Savin**  for  a  square  hole  with  rounded  corners 
and  somewhat  concave  sides  is  -1/6.  A  modified  version  with  fairly 
straight  sides  and  rounded  corners  can  be  obtained  by  taking  m^  ■  -1/8. 
For  comparison,  the  values  of  the  inegral,  I,  for  this  case  were  found 
to  be  0.3627,  0.3958,  and  0.4532  for  m^  =  -0.05,  -1/8,  and  -1/6, 
respectively.  Figure  1  compares  the  shape  of  this  optimized  hole 
(m^  »  -0.05)  with  the  shapes  of  holes  obtained  by  setting  m^  =  -1/8  and 
-1/6.  Figure  2  compares  the  boundary  stress  distribution  for  these  holes 
plotted  as  a  function  of  6.  The  highest  stress  is  smallest  in  the  case 
of  optimized  shape.  The  boundary  stresses  for  the  optimized  shape 
discussed  by  Durelli  and  Rajaiah^  are  also  included;  however,  a  complete 
distribution  as  a  function  of  8  could  not  be  plotted  as  the  relationship 
between  8  and  the  polar  angle  for  this  shape  is  not  known. 

Case  iii: 

The  third  case  (m^  =  -0.01,  p  =  2,  q  =  1)  represents  only  a  slight 
perturbation  on  a  circular  hole;  nevertheless,  for  this  stress  field,  it 
is  the  optimum  shape  for  a  hole  whose  boundary  is  constrained  by  a 
function  given  by  Equation  (13).  The  value  of  the  integral  for  this  case 
is. 1.3723. 

The  calculations  which  have  been  performed  for  Equation  (13)  were 
repeated  for  the  elliptical  shape,  i.e., 

z  =  C  +  —  *  lm1l  <  1  (18) 

The  value  of  the  integral  in  this  case  was  found  to  be 

(2Sm  +  D2) 

I  =  2S2  +  - =- -  (19) 

1-m^ 

which  for  a  minimum  value  requires 


It  is  interesting  to  note  that  this  value  of  m  coincides  with  that 
2  1 

of  Bjorkman  and  Richards  for  their  harmonic  hole.  The  value  of  the 
integral  in  this  case  is  simply  2S2,  i.e.,  1.125  when  p  =  2  and  q  =  l  and 
m1  =  1/3. 
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A  comparison  of  this  elliptical  case  with  case  (iii)  of  the  square¬ 
like  hole  (Equation  (13))  shows  that,  of  all  the  shapes  considered,  the 
elliptical  hole  produces  the  smallest  stress  concentrations  under  a 

2 

biaxial  stress  field.  On  the  basis  of  the  work  of  Bjorkman  and  Richards, 
this  conclusion  can  be  generalized  to  include  any  hole  shape.  However, 
if  the  hole  shape  is  constrained  to  a  shape  given  by  Equation  (13) ,  then 
case  (iii)  represents  the  optimum  situation.  Thus  whenever  m^  is  present 
in  the  mapping  function,  an  optimization  procedure  (without  additional 
constraints)  as  presented  here  leads  to  the  trivial  solution  of  all  m^ 
except  m^  equal  to  zero,  when  m^  is  given  by  Equation  (20).  Other  non¬ 
trivial  situations  could  be  investigated  by  imposing  additional  constraints 
such  as  a  priori  assignment  of  numerical  values  to  other  m^.  It  should 
be  emphasized,  however,  that  the  elliptical  optimum  shape  holds  true  only 
in  the  case  of  biaxial  stress  fields. 

Case  iv: 

The  fourth  case  (m^  *-.09,  p=l,  q  =  -1)  represents  pure  shear 
loading,  which  does  not  permit  an  optimum  elliptical  shape.  However,  the 
present  constrained  optimization  procedure  leads  m3  to  -0.09,  at  which 
point  I  =  0.9077.  This  optimized  hole  is  also  square-like  (double  barrel 
shape)  with  rounded  corners.  For  comparison,  the  circular  hole  in  this 
stress  field  will  produce  an  I  =  1.0.  The  maximum  stress  concentration 
is  3.07  at  8  =  23.6°  (8=0  corresponds  to  one  of  the  applied  stress 
directions)  which  compares  favorably  with  a  stress  concentration  of  4  for 
a  circular  hole  in  a  large  plate  under  pure  shear. 

EXAMPLE  II 

The  cases  discussed  so  far  dealt  with  the  mapping  function  containing 
only  one  term.  Inclusion  of  an  additional  term  in  the  mapping  function 
leads  to  considerable  complexity  in  the  algebraic  manipulations.  The 
mapping  function  considered  was 

m  m_ 

z  -  5  +  +  -f  (21) 

(-3  57 

In  this  case,  the  symbolic  expression  for  I,  the  stress  integral,  and  its 


derivatives  were  obtained;  however,  these  were  so  unwieldy  that  a  direct 
numerical  approach  was  followed.  It  was  found  that,  in  the  case  of  a 
uniaxial  stress  field  (p  ■  1,  q  =  0) ,  the  minimum  value  of  I  *  0.3623 
occurred  at  m^  =  -0.05  and  m^  *  0.0035.  It  is  interesting  to  note  that 
this  case  implies  an  improvement  over  the  shape  obtained  by  case  (ii)  of 
the  last  example.  The  actual  difference  between  the  two  cases  is, 
however,  very  small. 

CONCLUSIONS 

A  stress  integral  has  been  developed  which  can  be  effectively  used 
to  optimize  a  number  of  hole  shapes  in  structures  consisting  of  large 
plates.  For  example,  it  has  been  possible  to  obtain  a  stress  concentration 
factor  of  2.47  for  a  square-like  (double  barrel  shape)  hole  with  rounded 
corners.  This  value  is  about  3%  lower  than  the  optimized  quasi  square 
shape  and  about  5%  lower  than  the  double  barrel  shape  described  by 
Durelli  and  Rajaiah.^  Both  shapes  occurred  in  uniaxially  loaded  plates. 

In  the  case  of  pure  shear  loading  of  large  plates,  another  square¬ 
like  (double  barrel  shape)  hole  shape  has  been  found  which  reduces  the 
maximum  stress  concentration  to  3.07  from  a  value  of  4  for  a  circular 
hole. 
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APPENDIX  -  CONSTANTS  FOR  EQUATIONS  (A)  AND  (5) 


a^  ■  -  S 

a,.  *  ((15  m|  +  ((-9  -  5)m^  -  6  +  6  01^)111^  +  (3  ra^  m| 

+  (-  6  nij  +  nj  +  3)m,.  -  2  +  3  m|  +  +  (m^  -  l)m^)S  +  3  D  m^ 

-  D  m_)/(15  +  (9  m  -  5)m?  +  (-3  m  m,  -  m?  -  3)m  -  3  mf 

-  +  1) 

=  ((15  +  (-  10  m,.  +  9  +  5  “3)®^  +  (*  10  -  3  m3 

+  6  +  (m2  -  3)mj  -  2  m^)m^  -  3  +  (m^  -  2  m^)m^  +  m2 

-  ra^S  -  D  -  D  m^)/(15  +  (9  m^  -  5)m2  +  (-3  m^  -  m^ 

-  3)m^  -  3  aij  -  iij  -  +  1) 

=  ((15  m|  +  (30  +  9  m^-  5  +  ((-  3  -  10)m5 

-  6  -  ra|  +  3  m^)m^  -  3  mj:  +  (-  6  -  m^m^  “  m3 

-  m1)S  +  3  D  m?  -  D)/(15  m|  +  (9  n»3  -  5)m2  +  (-  3  m5  -  mj  -  3)m7 

-  3  m2  -  -  in^  +  1) 
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